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Abstract 

Some new developments in constrained Lax integrable systems and their applications 
to physics are reviewed. After summarizing the tau function construction of the KP 
hierarchy and the basic concepts of the symmetry of nonlinear equations, more re- 
cent ideas dealing with constrained KP models are described. A unifying approach to 
constrained KP hierarchy based on graded SL(r + n, n) algebra is presented and equiv- 
alence formulas are obtained for various pseudo-differential Lax operators appearing in 
this context. 

It is then shown how the Toda lattice structure emerges from constrained KP models 
via canonical Darboux-Backlund transformations. These transformations enable us to 
find simple Wronskian solutions for the underlying tau-functions. 

We also establish a relation between two-matrix models and constrained Toda lat- 
tice systems and derive from this relation expressions for the corresponding partition 
function. 
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1 Introduction 

1.1 Abbreviations; Diagram of Lectures 
Abbrevations used in text: 

KP hierarchy (equation) = Kadomtsev-Petviashvili hierarchy (equation) 

cKP hierachy = constrained KP hierarchy 

DB transfromation = Darboux-Backlund transfromation 

KdV hierarchy (equation) = Korteveg-de Vries hierarchy (equation) 

BA function= Baker- Akhiezer function 

AKNS hierarchy = Ablowitz-Kaup-Newell-Segur hierarchy 

AKS approach = Adler-Kostant-Symes approach 

GNLS equation= generalized non-linear Schrodinger equation 



Fig. 1. Structure of exposition 



KP Hierarchy 



symmetry 



constraint 



sl(r + n, n) cKP 



sl(r) KdV 



2-Matrix Model 



string 



equation 



DB 



transformation 
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1.2 Content of Lectures 



In these lectures we study constrained KP hierarchy, its connection to the discrete Toda 
models and application via Toda models to the multi-matrix approach to the two-dimensional 
gravity. 

A presentation follows the diagram shown on the previous page in the counter-clockwise 
direction. In Section 2 we collect some standard facts on KP hierarchy which are used in 
several places in these lectures. A link between the BA function and r function is estab- 
lished and it is shown how to obtain the conserved current densities of the KP hierarchy 
from these concepts. A complementary presentation of KP hierachy centered around multi- 
Hamiltonian/Poisson structures was given at the previous Swieca school ||. 

The notion of constrained KP hierarchies is introduced in Section 3. Since this notion 
is significant for our lectures we introduce a variety of reduction schemes and elaborate on 
their mutual equivalency. 

The picture which emerges from this discussion is most simply presented in the language 
of pseudo-differential operators. Two basic equivalent approaches look especially attractive 
and deserve to be mentioned here in the introduction. 

One construction introduces the Lax operator of the constrained KP (cKP ) hierarchy as 
a ratio 

L m n = —r^ n < m - 1 (1.1) 



J m,n 

of two purely differential operators 



LM = (D + v m ) (D + v m ^) ■■■(D + v 1 ) ; L<"> = (D + v n ) (D + v n ^) ■■■(D + v 1 ) 

(1.2) 

with coefficients V{,Vi subject to a constraint: 

m n 

E^-E^ = o (i-3) 

3=1 1=1 

Imposing condition (|1.3|) is equivalent to requiring tracelessness of the underlying graded 
SL(m,n) algebra. For this reason we refer to this class of constrained KP (cKP ) models 
as SL(m,n) cKP hierarchy. This construction is useful to describe a bi-Poisson structure of 
the cKP hierarchy and in fact involves variables which almost abelianize the second Poisson 
structure. There is another formulation of the same class of models which is based on an 
alternative expression for the Lax operator given in ( |1 . 1| ) : 

r—2 n 

L m=r+n<n = D r + £ uiD 1 + J2 QiD-^i (1.4) 

1=0 i=l 



A connection between variables from ( |1.1|) and (1.4) takes a form of complicated Miura-like 



transformations ( see Section 3 and [T^]), but for \l/j's the link is very simple since \l/j's turn 
out to be elements of the n-dimensional kernel of L^. Furthermore abelianize the 

first bracket structure. Their property of being eigenfunctions of L r+n ^ n makes the expression 
in ( |1.4| ) a proper framework for discussing the Lax equation in the setting of the SL(m,n) 
cKP hierarchy. 
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To go further and find r-functions corresponding to the SL(m, n) cKP hierarchy we need 
a notion of the Darboux-Backlund (DB) transformations which we introduce in Section 4. 
With this notion we are able to derive simple Wronskian expressions for the r-functions. 
We also reveal a connection between on one side (constrained) Toda discrete models and on 
another cKP models endowed with DB symmetry structures. The picture which emerges is of 
the continuous integrable model with canonical symmetry mimicking the lattice shift of the 
corresponding Toda lattice. This scenario has been appearing in various forms in literature 

The observations of Section 4 provide a right scene for establishing the link to two- 
matrix model. The two-matrix model is rewritten in Section 5 as a linear discrete system 
with additional constraint; a string equation. Using the string equation we are able to cast 
the lattice equations into the form of a constrained Toda equation which can be represented 
via construction of Section 4 to the class of cKP models. 

We have included two appendices on Schur polynomials and Wronskians, which may 
clarify few technical points. 



1.3 Acknowledgements 

These lectures have grown out of work I have done over the past few years with L.A. Ferreira, 
J.F. Gomes, E. Nissimov, S. Pacheva and A.H. Zimerman. I thank them for sharing their 
insights with me. 

It is a pleasure to thank organizers of the School for their warm hospitality and giving 
me opportunity to present these lectures. 
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2 Pseudo-differential Operators, BA Functions and Hi- 
rota Equations 

2.1 Lax Equation and Eigenfunction for the Lax Operator 

We begin with the Sato theory of the KP hierarchy. Let {tj} denote a set of independent 
variables with t\ = x. The formulation of the KP hierarchy is based on the Lax equations 

8L 

— = [B n ,L] n = l,2,... (2.1) 
describing isospectral deformations of the pseudo-differential operator: 

oo 

L = «9 + 5> i «r t - 1 (2.2) 

Here u n are the functions of {tj}, and B n = (L n ) + is the truncation to the differential part of 
L n . A negative power part of the differential operator L n will be written as (£")-. The Lax 
equations ( |2.1| ) have the zero-curvature representation taking the Zakharov-Shabat form: 

^-^ + [B n ,B m }=0, n,m=l,2,... (2.3) 

It is often convenient to view fl2.2|) and (2.3) as integrability conditions of the linear 

system: 

Lip(t,X) = \ip(t,\) (2.4) 

^ = iW(«,A) (2.5) 
dt n 

In connection with the above linear eigenvalue problem one introduces two class of functions. 
Definition. A function $ is called eigenfunction for the Lax operator L satisfying Sato's 
flow eguation (\2. 1\ ) if its flows are given by expression: 

<9$ 



dt. 



(L m ),$ (2.6) 



for the infinite many times t m . 

A special role is played by an eigenfunction function, which also satisfies the spectral 
equation (|2.4j ). 



Definition. A function ip(t, A) satisfying relations ( \2.4 ) and (2.5) is called a Baker- Akhiezer 
function. 

2.2 Conservation Laws of KP Hierarchy, and Tau-Function Rep- 
resentation of the Baker- Akhiezer Function 



Note first that (|2.2| ) can be inverted providing an expansion of D in powers of Lax operator: 

oo 

£ = ^ + E^ (1) ^ (2-7) 
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Similar inversion relation can also be written for the differential operator B m 

oo 

B m = L m + Y.°t n) L- i (2.8) 
i=i 

Expand now lnip(t, A) as 

oo oo 

A) = £U" + ( 2 - 9 ) 

n=l i=l 

Because of the linear equation ( |2.4j ) for L, one finds 

oo 

X(X) = m*, A) = A + E ^ (1) A- 1 (2.10) 

i=l 

where we have defined for convenience an auxiliary quantity x — din if). Hence by comparing 
(|3) and (U) we get ^ = af } . 

The conservation laws for the KP hierarchy in form we will discuss here follow from the 
recurrence relation for the differential operators B m . From Q2.8|) it follows namely 

(oo \ oo 

d-J2 U m + E ^ m+1) L- 1 
i=l / i=l 

Coo \ / oo \ oo 

D-Y, "PL-* )[B m - £ + £ af +1) ZT< (2.11) 

i=i / \ j=i / i=i 

m 

= DB m -Y vfBm-j - 4 m) + O (D- 1 ) 

3=1 

Since the left hand side is a pure differential operator it is obvious that terms with D~ k with 
k > cancel and consequently the only non-zero contribution comes from the first three 
terms on the right hand side. The following identity: 

[D - x(A)] B{\) — e(A) (2.12) 

with 

oo oo 

e(A) = -l+J2 A- m -Vi m) ; B{\) = E X~ m ~ 1 B m (2.13) 

m=l m=0 

presents a compact way of expressing the recurrence relation of equation (J2.ll). We now 



follow reference (see also |59], [yj) and apply (|2.12|) on the BA function ip{z) with 



expansion parameter being z instead of A: 



or 



E A-™- 1 ^^) - X (A) E A- m - 1 5 m ^(z) = e(A)V(z) (2.14) 

m=0 m=0 



E A— 1 (dB m ij;(z)) ■ ^\z) - X (A) E A^fWC*) ■ r 1 ^) = 6(A) (2.15) 

m=0 m=0 
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Using that B m ijj(z) ■ ip 1 (z) = z m + aiz 1 + ... we obtain 



e (A) = ]P X~ m ~ l d (B m ^{z) ■ - X{Z) _ * (A) (2.16) 



m=0 

oo 



z-X 



m=0 m=0 



which upon taking the limit z — > A yields 

<9%(A) 



, 6(A) = 9 £ A"™' 1 (5 m ^(A) ■ r\X)) (2.17) 

^ m=0 



Hence modes of the left hand side of ( 2.17 ): a® — la\ 1 ^ , 1 = 1,... are conserved densities of 
the KP hierarchy. From 

d {B m i> ■ r l ) = d(^- r 1 ) = (50 • r 1 ) = §- (2.18) 

we find that also da[ /dt m define conserved densities of the KP hierarchy and moreover 
comparing (|2.17|) and (|2.18|) we get a relation between these conserved quantities: 



E uu m~j (m) (i) / n , n N 

"op = CT 1 - mC7 m ( 2 -!9) 

i=i OT i 



From an obvious identity 



we see that we can write a[ m ^ = df/dt m with some arbitrary function /, which we will 
choose to write as / = — c^lnr. Hence in this new notation er^ = —d x d m lnr and clearly 
dn^i* 1 ^ — d m a[ . Furthermore we can now rewrite ( |2.19| ) as 

ft ftn^ 

-^lnr = -ma«-£^ (2-21) 



dt m k=1 dtk 

Remark. From ( p.7|) cr^ = — Resg(L) and therefore «o = d^lnr. Recall that for the 
pseudo-differental operator A = aiD" 1 + ... we have Res^A = a±. 

On basis of identity ( A.3f) (or ( |A.6|) ) for the Schur polynomials p n shown in Appendix A 



one verifies that ( |2.21| ) has a solution given by 

^^^(-^lnr ; n>l (2.22) 

where 

3=1 

Mi 2dt 2 ' Sdt 3 
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9=(4r,7^r,^r,-) (2-23) 



Comparing with and recalling that dipi = we see that (up to an integration con- 
stant) 

oo oo 

lnV(M) = X;*nA n + E^ _i Pi(-5)lnr (2.24) 

n=l i=l 

which yields the main result of this subsection: 

V(*,A) = e^ A > V Wtl) ^ =^)V ^ ^ — A"" (2.25) 

r W n=0 T W 

where 

00 

£(t,A) = 2*nA" (2.26) 



n=l 



2.3 Dressing Operator 



It is possible to reproduce the Lax operator (2.2) through the dressing formula 



L = WDW- 1 (2.27) 

where the dressing operator W is the pseudo-differential operator: 

00 

W = 1 + 2 (2.28) 

satisfying the Sato equations: 

dW 

= B n W - WD n = - (L n )^ W (2.29) 

ot n 

LW = WD. (2.30) 

From the dressing formula it follows immediately that the BA function can be rewritten as 

00 

jj,(t, A) = Wexp(J2 ^A n ) = w(t, A) e 5(t ' A) (2.31) 

?1=1 

where 

*»)=^ (2.32) 

as can be found from fl2.25|) . 

In a similar way an adjoint BA function i[>*(t, A) can be constructed as follows: 

i/j*(t,\) = (W l Ye~^ x) =w*(t,X)e-^ (2.33) 

= iMM (2 . 34) 
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and defines a linear system: 



L*il>*(t,\) = \i/>*(t,\) (2.35) 

dr(t,\) 



dt n 



-BW\t,\) (2.36) 



where L* is adjoint of L and B* is a differential part of (L*) n . Recall that the adjoint 
operator A* of A = aD k + ... is given by A* = (—l) k D k a + 



2.4 Hirota Equations 

As seen from (p.31|) and ( |2.25| ) we can rewrite the dressing operator as 



leading to 



L n = WD n W~ 1 = £ Pjt^pn-^Pi^m (2 . 38) 
Let us insert the above expressions into the Sato equation ( [2.29|) . Noticing that 

_ m _ = _"+'g + ' tMir. + . . . (2.39) 

i,j>o T W 



and taking the residue on both sides of fl2.29|) we get: 



or 

'iD 1 D n -p n+ i(5))r-r = (2.41) 
where we used Hirota's operators defined by 



D T a ■ b = + s Mtj ~ *i)l.,=o (2.42) 



The first non-trivial eqution obtained from ( |2.41| ) for n = 3 ( by dropping odd polynomials 
in D, which do not make any non-zero contribution) is the KP equation: 

(p\ + 3Dl - 4 J D 1J D 3 ) r • r = (2.43) 

which in components takes the usual form: 

d_ fduo _ 1 Ohio , ^o^\_3^o_ n (0 AA ^ 

dx { dU 4 dx* dx Adt 2 ~ 1 > 



with uq = ResgL = d% In r. 

Remark. Eq. (|2.41|) is contained in the Hirota differential equations for the KP hierarchy 
taking the following expression: 



52pj(-2y)p j+1 (D t ) exp < E^A ; 
j=o u=i 



r(t) ■ r(t) = 



(2.45) 



with y = (|/i, y 2 , . . .) being an extra multi- variable. 

Indeed, it is easy to see that (|2.41|) can be obtained from fl2.45|) as a coefficient of the term 
linear in y n . It can be shown that the Hirota differential equations ( |2.45[ ) can be obtained 
from the bilinear identity for the r functions: 



T t 



1 1 



(2.46) 



by expanding arguments of both r-functions in ( |2.46| ) [p4], 59|. Hence ( |2.46| ) carries a complete 
information about the KP hierarchy. The KP formalism based on ( |2.46|) instead on the 
pseudo-differential Lax operators can be analyzed in terms of the free fermions Fock space. 



These beautiful results are beyond the current discussion (see for details p4|, [39]). 



2.5 Hamiltonian Densities and the r-Function 



Let us now consider the product s = ipip* of the BA function and its adjoint. We will now 
show that s x generates infinitely many commuting KP flows. First as in note from Q2.32|) 
and pi) that 



Tjti-l/iX^Tjti + l/iXi) 

r 2 (U) 

1 ~ p k (p) 



r 2 (U) 



exp [Y^ 







\ k\ k de k 



t (t + e) r (t - e) 



e=0 



r • r 



E 



DxD^T ■ T 



2r* (t) ^ ^ 



(2.47) 



where in the last equation we have used the Hirota equation fl2.41| ). This equation suggests 
to consider the Laurent coefficients of s defined through 



n=0 



DiDn-iT ■ r 
2r 2 



(2.48) 



We now recall the following technical lemma |2(J : 

Lemma. Let P, Q be pseudo- differential operators. We have 



Res> 



{Pe xX ) (Q 



* e -x\ 



Res d (PQ) 



as follows from a direct verification. 

Taking P = WD n and Q = W~ x we obtain from the Lemma that 



(2.49) 



s n+1 = Res A [X n W\ = Res a (L n ) 



(2.50) 
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From (1217| ) we find 



1 



DxD n -]T ■ t = d x d n -i lnr 



(2.51) 



2r 2 {t) 



Since mq = <9 2 In r we can rewrite d x s n as 



9m, 



o 



(2.52) 



dt 



n— 1 



Note that in notation of the previous subsections the above equation shows that the conserved 
densities of the KP hierarchy a[ m ^ are equal to Hamiltonian densities 7i m = Resg (L m ) of the 
KP hierarchy, which generate commuting flows through relation d x 7i n = (d/dt n )uo- Accord- 
ingly the KP hierarchy has an infinite set of commuting symmetries associated with infinite 
set of the conserved Hamiltonians. This translates into commutativity of Hamiltonians on 
the level of the Poisson brackets indicating integrability of the system. These symmetries are 
called isospectral symmetries as they preserve the spectrum of the underlying linear problem. 
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3 Symmetry Reduction of the Integrable Systems. Five 
Constructions of the cKP models. 



3.1 Symmetry of Nonlinear Evolution Equation 

We will discuss notion of symmetry of partial differential equatiom on the basis of an evolu- 
tion equation: 

with some nonlinear operator K. An evolution equation |^ = G[u(t)] is called a symmetry 
of fl3.1| ) if the Frechet derivative of K at the point u in the direction of G 

K'[u{t)](G) = ^-K[x,u(t) + eG,u' + eG x ,..]\ e=0 (3.2) 

o. = !+gHr"* +1 » (3.3) 

satisfies 

K'[u{t)] (G) = + G'[u(t)] (K) (3.4) 

Hence a condition that G is a symmetry of ( |3.1|) is equivalent to having v = G[w(t)] satisfy 
the linearized version of ( |3.1|) with a background u(t) i.e. dv/dt = K'[u(t)] (v). In other 
words u + ev satisfies ( |3.1|) for all solutions u of ( p.l|) for an arbitrarily small e. Therefore the 
linearized version of the nonlinear evolution equation contains all the pertinent information 
about its symmetries. 

We can also rephrase condition (|3.4|) as commutativity [ D t , D a ] = of two differentia- 
tions associated with t and a 

°' - l + A l + ^ + ... (3.5) 

d d d 
D a = — + G— + G X — + ... (3.6) 
oa ou au x 

Examples: 

Consider the KdV equation: u t = 6uu f — u'" in terms of uj' = u: 

u t = K[u] = 3u' 2 - J" (3.7) 

Condition (|3.4|) gives 

dv/dt = K'[u{t)\ (v) = 6u'v' - v"' (3.8) 
Consider Galilean symmetry of KdV: 

uj — ► uj = uj(x + 3te) H — xe (3.9) 
which can be cast in the form of "evolution" equation: 

v — G — ^T-\ € =o = 3tuj' + -x (3.10) 
de 2 
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Taking into account that v x = 3tuj" + 1, v xxx = 3tu IV and also dv/dt = 3uj' + 3t [3uj 12 — uj'") 
we can verify that the Galilean transformation is a "non-isospectral" symmetry of the KdV 
equation. 

For the KP eq. (|2.44 ) its linearized version takes the form: 
3.2 Symmetry Reduction and the cKP Hierarchy 

We will first show that the quantity <9 X ($\I/) is a symmetry of the KP equation with $ and 
\l/ are (adjoint ) eigenfunctions of L i.e. they satisfy (|2.5| ) and ( |2.36|) : 

-$ = L™$ ; — $ = -L m ;^ (3.12) 



dt + dt 

where L m * is the adjoint operator of L m . 

Let us particularly stress that the above eigenfunctions do not need to be Baker- Akhiezer 
eigenfunctions of L and d x (Qty) in general differs from quantity s x = d x ipip* , which, as 
we have seen in the previous section, generates an infinite set of commuting symmetries 
associated with infinite set of the conserved Hamiltonians of the KP hierarchy. 

For n = 2 equation ( |3.12| ) gives: 

<9$ <9 2 $ d^ d 2 ^ 

wr™ +2 ^ ; wr-^- 2 ^ (3 - 13) 

while for n = 3 we obtain from ( |3.12| ): 

d^ d 3 V o <9# 3 /T 3 f 1 duo\ /01/n 

^ = ^ +3m °^ + 2^-2^ w 2 r (3 - 14) 

Using equations (|3.13|) and (|3.14|) one can show as in |l8j that S = $ 1 I / satisfies the 
evolution equation 

ds id*s as a^sPs* 

df 3 = ^ + 3u °d^ + f ~W ( } 

and therefore the quantity S x = d x S satisfies a linearized version of the KP equation (|3.11| ). 
Hence according to the above definition S x is a symmetry of the KP equation. 

Recall now that in components the Lax equation (|2.1|) takes the following form of partial 
differential equations 

d r U n = F n>r (u Q , Mi, ... , Mr+n-l) 71 = 0, . . . (3.16) 

with differential polynomials F nr in u Q ,ui, . . . , w r+n _i. In particular we have a hierarchy of 
flow equations: 

d r uo = F r (u ,u 1 ,...,u r -i) (3-17) 
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Equations ( |3.17| ) can be viewed as symmetries of the KP equation according to our above 
definition. 

We will first introduce a constrained KP (cKP ) hierarchy by imposing a constraint on 
symmetry represented by the evolution equations in fl3.17| ). 

We will discuss three (very related) approaches to define the cKP hierarchy and show 
their equivalence. 

Symmetry constraint, cKP hierarchy, Definition I 

This approach defines the constrained KP hierarchy by imposing the symmetry constraint 
on flows in ( |3.17| ) in terms of the eigenfunction and its adjoint entering the linear problems 
(PD and (ET5ED : 

F r {u , it r _i) = (3.18) 

Hence for each r = 1, 2, ... we define the different cKP hierarchies by imposing equality of 
the flows/symmetries F r to the symmetry d x (<&ty) of the KP equation. The constraint allows 
us to eliminate the Lax coefficients Uk> r -i in terms of eigenfunctions $ and \& and «/. <r _i. 

The next two approaches will allow us to write down an explicit Lax representation of 
the cKP hierarchy. 

Symmetry constraint, cKP hierarchy, Definition II 

Consider the Lax operator Q = D + J2tLo UiD^ 1 ^ 1 satisfying the flow equations as in ([2.1Q 
i.e. 

^- = [B k ,Q] ; B k = (Q k ) + k = l,2,... (3.19) 

Let furthermore (like in fl3.12|) ) $ and ^ be, respectively, eigenfunction and adjoint eigen- 
function of Q, meaning that for k = 1,2,... it holds 

= B k § ; — = -B* k V (3.20) 
dt k dt k k y 1 

where B^ is an adjoint operator of B k . We now impose the constraint on the KP hierarchy 
by requiring that the Lax operator Q satisfies a condition [ pO|j : 

Q r = B r + ^ZT 1 ^ ; B r = (Q r ) + (3.21) 

for r being a fixed positive integer. 

Let us first verify that Definition II implies the symmetry constraint of Definition I. Let 
us namely impose constraint ( [3.2 1|) and notice that from ( |3.19|) we get 

du /dt r = Res 9 (dQ/dt r ) = Res d [B r , Q] = Res d [Q r - ^ZT 1 ^ , Q] (3.22) 
= Resap^Zr 1 ^, Q] = d x ($V) 

This clearly shows that the Definition II realizes directly in the Lax setting the symmetry 
constraint of Definition I. 

Lemma. The time evolution of the pseudo- differential operator ^D^ 1 ^/ is qiven by: 

d 

—QD^m = [B k , ^ZT 1 ^]- (3.23) 
ot k 
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The proof is a consequence of following technical observations based on ( |3.2(J| ) 



(B k $D- 1 y)_ = (B k <t>) D-^ = gfD- 1 * 



and 



^d~HBi 



dt k 



(3.24) 



(3.25) 



What we learn from the above lemma is that if we define the Lax operator L such 
that its purely pseudo-differential part is L_ = $>D~ 1 ^ then L_ satisfies automatically the 
KP flow equations: dL^/dtk = [Bk , Note that in the last equation we used that 

[Bk , L\- = [Bk , L_]_. These observations lead us to the next definition of cKP hierarchy. 

Symmetry constraint, cKP hierarchy, Definition III 

Here we work with Lax operator of the cKP hierarchy defined as 

r-2 



dt k [ V 



,L] 



(3.26) 
(3.27) 



We first address an issue of equivalence between Definition II and Definition III, which is 
easy to establish when making a connection Q r = L. It is then easy to see that the two flows 



given below are equivalent (B k = (Q k ) = (L^ 1 



dQ/dt k 



Q k ), , Q 



dL/dt k 



L 



L k/, 



(3.28) 



We have already seen in ( [3.23D that the pseudo-differential part of L from ( |3.26|) evolves 
according to the Sato's formalism. Let us now investigate evolution of the differential part 
B r = (Q r )+ = (L) + of L: 



dB r /dtk = [B k , L} + = [B k , B r ] + [B k , QD-^l 
which after using Zakharov-Shabat equation (|2.3| ) becomes: 

dBk/dt^^k^D- 1 ^ 



(3.29) 



(3.30) 



valid for all k = 1,2,... . Equation ( |3.30| ) describes how the cKP constraint condition is 
being imposed on the flows of KP coefficients in Bk. Let us study consequences of (|3.30|) in 
some simple cases. The first nontrivial case occurs for k = 2 with B2 = {Q 2 ) + = D 2 + 2uq = 
D 2 + U . From flOGD it follows: 



2du /dt r = [D 2 , $ J D" 1 ^]+ = 2d x (m) 



(3.31) 



which reproduces the flow equation being basis for the Definition I of the cKP hierarchy, 
establishing henceforth equivalence between Definition I and Definition III. 
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For k = 3 we have B 3 = (Q 3 
( p.30| ) it follows this time that: 



= D 3 + 3u D + 3m + 3d x u = D 3 + U D + U x . From 
d(3u D + 3ui + 3u' )/dt r = [D 3 , SZr 1 ^ ] + = 3$"^ + 3 W + 3 ($*)' D (3.32) 



or 



duxjdtr = - ($*')' (3.33) 

Consider the Lax operator L having the form as in (l3~26j ) and satisfying ( gg7p . One can 
ask whether $ (\I>) are automatically the (adjoint) eigenfunctions. The following Lemma 
answers this question affirmatively 

Lemma. Let L = L + + $.D _1 \I/, with r being the order of L + . Then the Lax equations 
of motion -^-L = [B^ , L\ imply: 



wiiere 5 fc = (L fc / r )+ , 5£ = (L fc / r ); 

Remark. Eq. (|3.34 ) is equivalent to the infinite set of equations: 

(^*-w*)„)^*+*4(^*+(W 

An obvious solution of ( |3.35|) seems to be: 



(3.34) 



Z = 0,1,2, 



(3.35) 



(3.36) 



where c r is x- independent. Thus, up to a x-independent phase transformation $ —>■ e Cr $ and 
^ — > e _Cr \l/, which does not change the Lax operator, $ fi?) are (adjoint) eigenfunctions. 



Examples: 

r=l. The simplest case is when r = 1 for which we have: 

F x = d x u = d x ®^ -> m = + const 
Inserting this into ( |3.13| ) (with zero integration constant) we get 

d 2 m 



= h 2$ $ ; 



- 2$^ 2 



(3.37) 



(3.38) 



in which we recognize the Nonlinear Schrodinger (NLS) equations being the first (non-trivial) 
flow of the AKNS hierarchy. Eqs. ( |3.38| ) are also results of the Sato equation ( p.lQ for two- 
boson hierarchy |4l| , |5| defined by the Lax operator Q = L = D + QD^ 1 ^. Expressions for 
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the coefficients Uk>o of the standard expansion for Q are easy to find in terms of the Faa di 
Bruno polynomials || [| . 

r=2. Here F 2 = o 2 uq + 2d x u\ = d x which allows to express u n >i by uo, $, ^. 

Alternatively writing in the spirit of the Definition III the Lax operator as L = Q 2 = 
D 2 + 2u + ^D- 1 ^ we obtain from the Sato eq. (|3~28|) or (|3~T^ ): 

+ 2w $ 



d<5> 

dt 2 
duo 
dU 



dx 2 
d 2 ^ 
dx 2 



(3.39) 



which agrees with first flow equation of the so-called Yajima-Oikawa hierarchy [^T, [40 

r=3. Here F 3 = d 3 uo + 3d 2 u\ + 3d x u 2 + Quq3 x uq = d x ($^) which allows to express u n > 2 
by Uo, Ui, \&. Alternatively, writing in the spirit of the Definition 111 the Lax operator as 
L = Q 3 = D 3 + 3u D + 3ui + 3d x u + ^D' 1 ^ we obtain from the Sato eq. (ET2Sj) : 



9$ 

duo 
dt 2 
dui 



d 2 $ 

dx 2 

d 2 m 



dx 2 



+ 2w $ 



d 2 u + 2d x u 1 



-dim - -d' 3 u 



(3.40) 



2u d x u + -d x 



in which one recognize the t 2 flow of the so-called Melnikov system [fig] . The flow eqs. for 
t 3 = t r follow from (|3T31| ), ( ggg ) and ( |3~T2| ) (or (gig)): 



9$ 



-^r + 3W ^- 

ox 6 ox 



3«i + 



<9 3 ^ 

— r + 3u ^ 3«iW 



(3.41) 



In the next subsection we will learn how to view the above examples as special cases of 
SL(r + 1, 1) cKP hierarchy with r = 1, 2, 3 respectively. 

Remark. In the above we have kept for simplicity discussion restricted to a single pair 
$, \1/ of (adjoint) eigenf unctions. It is very simple to generalize all formulas to the case of n 
$ n , \l/ n (adjoint) eigenfunctions. Clearly in expressions ( |3.21| ) and ( |3.26 ) we need to make a 
substitution: $.D _1 \I/ — > J27=i ^iD~ l ^i- This is the structure, which will arise naturally in 
the next subsection. 
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3.3 cKP Hierarchy as SL(m, n)-type Lax Hierarchy and its Bi- 
Poisson Structure 

Let us introduce a ratio of two differential operators (with m = r + n in notation of the 
previous subsection). 



r = 



i» 



n < m — 1 



(3.42) 



where 



L (m) = D m + u m - 1 D m - 1 + ... + u l D + u ; = D n + u r ^ 1 D n ~ 1 + . . . + u x D + u (3.43) 



(«) 



} be a basis for the kernels of and 



Let {-ipi , i — 1, . . . , m} and \ipl ' , i — l, . . . , n 

L^ n \ respectively, i.e. we have L^ip^ = , % — 1, . . . , a — m or n. Alternatively, we can 
rewrite ( |3.43|) as 

= (£> + Um ) (£> + ...(£> + Vl ) ; = (£> + u n ) (£) + {}„_!) ■■■(D + V 1 ) 

(3.44) 

with 

^_ 1 [^ m) ,...,^ 1 Y 



Ui = d log ■ 



Wn = l 



and similar expressions for Vj. 

Taking the above into consideration we can rewrite (|3.42 ) as [ 12 , |B"2 



(3.45) 



£m,n = H {D + V 3 )\{{D+V l 



(3.46) 



1=1 



(see also |16], [25], [53| for discussion of the pseudo-differential operators of the constrained KP 
hierarchy). We will now study the Poisson structures of the hierarchy defined by the Lax 
operators of the form given in ( [3 .46 ). It is easy to prove the following Proposition |62|, [12]. 
Proposition. The Poisson bracket relation: 

{vi, v j } PB = 6 ij 5'(x-y) i,j = l,...m 

{v r ,v s } PB = -6 rs 5'(x-y) r, s = l,...,n (3.47) 
{Vi, vi} PB = 

is equivalent to 

{ (£"m,n\ X^ , (/^m,n| } PB A ( (^m,n-^0_|_ £m,nY ^- , m,n}^. Y C m n*j (3.48) 

where the subscript PB stands for the Poisson bracket as defined by ( 3.47|) . Moreover the 
Adler trace {C m>n \ X) defines an invariant, non-degenerate bilinear form: / cixRes (C m>n X) 
for the pseudo-differential operator C m>n and its dual X (see for instance || [l! 
We now introduce a Dirac constraint: 



j=i 1=1 



(3.49) 
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which is second-class due to: 

{# m , n , ^ m ,n}pB = (m - fl) 6'(x ~ Jj) 



(3.50) 



Condition (|3.49|) expresses tracelessness of the underlying graded SL(m,n) algebra. Cor- 
responding Poisson algebra of diagonal part of the graded SL(m,n) Kac- Moody algebra is 
obtained by the usual Dirac bracket calculation: 



{Vi, Vj} DB 
{v r , V S } DB 
{Vi , Vr}DB 



$ij - 



m 



- — )s'(x-y) i,j = l,.. 
— nj 

— — )$'(x-y) r,s = l, 



, m 



n 



(3.51) 



m — n 



5'{x-y) 



Proposition. The Dirac bracket ( ]3. 5\ ) takes the following form in the Lax operator rep- 
resentation: 

{ \£m,n\ X/ > \^m,n| Y/ 5 jjb \ \^wi»n| -^0 > \^m,n| ^/ } pb 

+ ! dxRes([£ m , n , X])9- 1 Res([£ m , n , Y]) 

in — a .1 

The proof follows from evaluation of the extra term of the relevant Dirac bracket: 

J { (£m,n | X) i ^ m.n } pp { ^ 'm.n i ^ 'm,n } p,B { ^ 'm.n > (£"m,n \ Y'j } (3.52) 

One easily verifies ( |3.52j ) using ( |3.50| ): 

{ (C m>n \ X) , * m , n (z) } PB = - - z) , £ m , n ]| X) = -Res ([X , C m>n ]) {z) (3.53) 
Formula (|3.52 ) contains as special cases n = m — 1 corresponding to the constrained KP 



hierarchy with L + = D and n = corresponding to the m-KdV hierarchy (see e.g. |27 



For the intermediary cases < n < m — 1 this formula presents a compact expression for 
the bracket structure of the SL(m,n) cKP hierarchy. 

The Lax operator from ( |3.42|) restricted to the constrained manifold defined by ( |3.49|) 
can be parametrized as 



n+l 



m— 1 \ n 



w = ^ u m , n=0 = n p+5)ii( fl +5H)p-E^ npHr 1 (3.54) 



l=m— 1 i=ra 

We can alternatively rewrite the expression ( |3.54j ) as 



z=i / ;=i 



n n 



m—n—2 

J2A l H(D + v i )- 1 + Yl A l+n+l D l + D m - n 



(3.55) 



i=i 



1=0 



with the second bracket structure automatically given by the formula ( |3.52| ). 
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For the special case m — n = r = 1 we have a canonical representation for variables (vi, vj) 
with 1 < i < m , 1 < I < m - 1 in L m , m _i = \~[) =m (D + v,-) \~\Z~i 1 (D + vi)' 1 in terms of 
the pairs (c r , e^J^^ , which are the "Darboux" canonical pairs for the second KP bracket 
satisfying {ci(x) , &i{y)} = —5ud x 5(x — y) for i, I = 1,2, ... ,m — 1. The representation is 

0: 



-ei_i 



to— 1 
m— 1 

- E c p 



/ = 1,2, 



"Cj + Uj 



, m 



(3.56) 
(3.57) 



Eq. ( |3.57|) includes the special cases of v\ = — Y^=i c p an d v m = —e\. One checks easily 
that in this representation the Poisson algebra of vi,Vi takes indeed the form of the bracket 
algebra of graded SL(m,m — 1) Kac-Moody algebra in a diagonal gauge: 



{vi, v j } = (5 lj - 1) 5'(x-y) i,j = l,...,m 

{v P ,v t } = - (8 p i + 1) S'(x -y) p, I = 1, . . . , m - 1 

{v l} vi} = -5'(x-y) 



(3.58) 



We can interpret ( p.54|) as a superdeterminant of the graded SL(m, n) matrix in a diago- 
nal gauge, which for KdV case n = becomes an ordinary determinant as in Fateev-Lukyanov 
[ 29| expression: 



J mS) 



D m + A m ^D 



m-2 



i=i 



m 

i=i 



(3.59) 



Example . As an example let us take m = 2, n = 1 in (|3.54|) . Then: 

L 2 ,i = (D + Cl + vx) (D - Cl ) (D + v,)- 1 = - [({?! + Cl )' + c x (v x + a)] (D + v^' 1 + D 



(3.60) 

as follows by direct calculation making use of an obvious identity (D + c) (D + v)^ 1 = 
(c-v) (D + v)- 1 + 1. 

One notices the absence of the term proportional to ]j m - n ~ x in ([3.55|) due to the SL(m, n) 
trace zero condition. Another feature is the presence of the constant term A n+X in the Lax 
operator L m n (for m — n> 1). This fact enables us to prove that SL(m, n)-cKP hierarchy is 
a bi- Poisson hierarchy. Consider namely L' m n = L m>n — A obtained by redefining the D° — 1 
term in the Lax operator by addition of the constant A. Clearly the value of the of bracket 
( |3.52j ) for the new Lax is 

{ ( L 'm,n X ) ' ( L m,n Y ) } DB = Tr ^ ^ {^'m,n X ) + L 'm,rX ~ ( X L 'm,n) + Y L 'm, 
1 



+ 



m — n 



dxRes ( 



±J m,n i 



d^Res ( 



r> y 

m.n ' 



X(L'\[X,Y] R ) (3.61) 



where we introduced an i?-commutator [X, Y]r = [X+, Y + ] — [X_, Y_] with subscripts ± de- 
noting projection on pure differential and pseudo-differential parts of the pseudo-differential 
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operators X, Y. Define next an i?-bracket {•, -}f as a bracket obtained by substituting 



^-commutator [X, Y]p> for the ordinary commutator |54], |58|, |10 



{ (L\ X) , (L\ Y) }f = — (L | [X, Y] R ) (3.62) 

Relation ( |3.61| ) shows that {•, -}db + A{-, -}f satisfies the Jacobi identity. We can state this 
result as: 

Proposition. SL(m, n)-cKP hierarchy is bi-Poisson with brackets {•, -}db an d {•, -}f defin- 
ing a compatible pair of brackets. 

This Proposition establishes the fundamental criterion for integrability of the SL(m,n)- 
cKP hierarchy. Clearly the argument holds also for the case m — n = 1, where one adds the 
constant A to zero representing the missing constant term. 

Let us now concentrate on the pseudo-differential part of L m n from ( |3.55|) . We can 
rewrite it as 

n n 

(L m , n )_ = Y, r iU D ~ 1( lr (3-63) 

1=1 i=l 

r^Ap-f* ; q n = el^ , gi = e/(*-*+0 , t = l,...,n-l (3.64) 
Let us define the quantity 

Q M = (-I) 4 "" fqj q^ [...[% {dx')^ 1 l<l<i<n (3.65) 



Then using that D = D l Qi ti D — Qx ti we obtain for quantity in eq. Q3.63|) 



(W)- = E r i 1] II D ' l( H + n D ~ X (-Qi,«-i9n) (3-66) 

where 

rf ' = n + riQ 1(i _i i = 2,...,n (3.67) 
The above process can be continued to yield an expression 

n 

(L m>n )_ =^$ i D- 1 ^ (3.68) 

i=i 

with 

i— 1 i— n i—n i—n 

= r i + E r ™ E ' ' ' E E Qn,i—8i-n-i—lQi—Si-n-i,i—8i-n-2—l 

n=l Sj_„_i=Sj_„_2+l s 2 =si+l si=l 

• • • (5i_ S2ii _ Sl _i(5i- Sl) i-i 1 < i < n (3.69) 

V n = q n , *i= {-l) n - l q n J q n ^ J ... J ^{dx')^ 1 < i < n - 1 (3.70) 

Note that the new variables coincide with elements in the kernel of L^ n \ It follows 
from the construction opposite to the one shown above and involving the following relation 
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(see for more details): 

n 



i=l 

n 



E4 W (D + B 



(3.71) 
(3.72) 



i=l 



where the new variables are : 



-4 



(n) 



5 



to 



(-i) M E* 

cLln 



W[* n ,...,* m ,*< 



, ^[*, vi/,., 



w[*n,...,* i+ i] 



(3.73) 
(3.74) 



From ( |3.72| ) we see that Vi = Bf 1 ' and from ( |3.74| ) it follows that ijj^ = 
Example . For n = 2 relations (|3.73|) - (|3.74f) become 



4 2) 
4 2) 



$^ + $2*2 , 4 2) = (^ln^- 1 ^ 
a x .ln^ 2 , b{ 2) =d x \n[y l (d x \n^ 1 y 



(3.75) 
(3.76) 



Equivalence between ( |3.72j ) and Q3.71| ) follows then by inspection. 

Proposition. $j, ^ are canonical-Darboux fields for the first bracket of the SL(n + 1, n) 
cKP hierarchy : 

{ $, ; , ^ } x = -<^(> -2/) , i,j = l,...,n (3.77) 

3.4 Affine sl(n + 1) Origin of SL(n+ l,n) cKP Hierarchy. 

In this subsection we will establish a connection between Generalized Non-linear Schroedinger 
GNLS matrix hierarchy for the hermitian symmetric space sl(n + 1) p3| and the constrained 
KP hierarchy. 

We first introduce ZS-AKNS scheme arising from the linear matrix problem for A £ Q 



with Q being a Lie algebra |32|, |49 



d tm V = B m y m = 2,3 



f> 9 ,l 

d= drr dx 



(3.78) 
(3.79) 



The compatibility condition for the linear problem ( |3.78| ) and ( p.79| ) leads to the Zakharov- 
Shabat (Z-S) integrability equations 



d m A - dB m + [A , B m ] = 
Let us use the following decomposition in (|3.78|) : 

A = XE + A with 



E 



d m = d u 



of. 



(3.80) 



(3.81) 
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where \i a is a fundamental weight and a a are simple roots of Q. The element E is used to 
decompose the Lie algebra Q as follows: 

Q = Ker(adE) © Im(adE) (3.82) 

where the Ker(ad-E) has the form JC' x u(l) and is spanned by the Cartan subalgebra of Q and 
step operators associated to roots not containing a a . Moreover Im(ad.E') is the orthogonal 
complement of Ker(ad-E). 

From now on we consider Q = sl(n + 1) with roots a = ctj + oti + \ + . . . + a j for some 
i,j = l,...,n,E= 2 ^ H ; ^ n is the n th fundamental weight and if a , a = 1, . . . , n are the 
generators of the Cartan subalgebra. This decomposition generates the symmetric space 
sl(n + l)/sl(n) x u(l) (see [f| for details). 

ZS-AKNS scheme becomes in this case the GNLS n (=sl(n + 1) GNLS) hierarchy. In 
matrix notation we have: 

/ 1 \ 
1 

1 



E 



1 



71 + 1 



V 



-n 



(3.83) 



+ r aE-(a a 



■+a n ) 



(3.84) 



The model is defined by A G Im(ad-E') and is parametrized by fields q a and r a , a — 1, . . . , n 
according to: 

A° = 

a=l 

which in the matrix form can be rewritten as 

/ • • • • • • q x \ 
••• q 2 

'■• : 



A 



(3.85) 



\ ri r 2 • • • r n / 
It can be shown that the flows as defined in (|3.78 ) and ( |3.81[) satisfy the recurrence relation: 



d m A° = Tld m ^A c 



K 



d — ad^o d ad ao ) ad 



A° 



(3.86) 



where we have defined a recursion operator 71 [[33|, (J. 

The connection to the cKP hierarchy is first established between linear systems defining 
both hierarchies. With (|3.83 ) and (p.85|) the linear problem from (|3.78 ) is explicitly given 
by: 

f d - A/(n + 1) • • • q 1 \ ( \ 

d-X/(n + l) ••• q 2 







r 2 



d-\/(n + l) q n 

r n d + n\/(n + 1) J 



'lf-2 



o 



(3.87) 
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Perform now the phase transformation: 
ipk — > = exp ' 



n 



Xdx) ij)k 



k = 1, 



n + 1 



(3.88) 



We now see that thanks to the special form of E in A = XE + A , ( |3.87| ) takes a simple and 
equivalent form: 

( d • • • qi \ { $i \ 
d ■■■ q 2 



i>2 







(3.89) 



9 q n ip n 

\n r 2 ■■■ r n d + A j \ ip n+1 J 

The linear problem ( |3.89| ) after elimination of ^ , k = 1, . . . , n takes a form of the scalar 
eigenvalue problem: 

n 

- [d - r k d~ 1 q k ]ip n+1 = Xi> n+1 (3.90) 
fc=i 

in terms of a single eigenfunction ^„ + i and the pseudo-differential operator 



L n = 9 - r k d 1 qk 

k=l 



(3.91) 



This formally relates GNLS n hierarchy to the SL(n + l,n) cKP hierarchy on basis of cor- 
respondence of the linear problems characterizing them. To fully establish their complete 
equivalence one can show that both models possess the same recurrence operators and there- 
fore all the flows are identical see || (|20| describes the the case n = 1). We first note that 
the successive flows ( |3.86| ) related by the recursion operator ( |3.86| ) are given by 



d n 



qi 



(i,l),(j,m) 



T 3 

qm 



(3.92) 



d, 



n-l 



( (-8 + r k d 1 q k ) + r { d 1 q j rid l r m + r m d 1 r i 

\ -qid^qj - qjd^qi (d - q k d~ l r k ) 5i m - qid~ l r m J ' ' \ q v 

From relation between the recursion matrix and two first bracket structure 1Z = i^-Pf 1 an d 
( fj.9'2j ) one finds an explicit expression for the second bracket to be : 

p _ ( tid^Vj + rjd~ l ri (d - £ fe nd^qk) S im - Tid' x q m \ , . 

2 V ( d -T,kqkd~ 1 r k )6ij -qid^rj qid' 1 q m + qmd^qi J 

in the same basis as in ( |3.92j ). 

The following Proposition proves the equivalence between sl(n + 1) GNLS hierarchy de- 
fined and the SL(n + l,n) cKP hierarchy introduced in the previous subsection. 
Proposition. Flows 

[«) + , Ln] (3.94) 



9t m L n 
dt m qi 



d-^r k d l q k 



k=l 



(3.95) 
(3.96) 



k=l 
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of SL(n+l,n) cKP hierarchy containing the Lax operator from ( 3.91 ) coincide with the flows 
produced by the recursion operator 1Z ( S.93J of the sl(n + 1) GNLS hierarchy 

Proof is given in || and is accomplished by showing that both hierarchies have identical 
recursion operators. Especially ( |3.96| ) yields for m = 2: 



% 
dt 2 
dr\ 
dt 2 



-d 2 x qi + 2q { Qb H 



6=1 



Efyn - 2ri Qb r b 

6=1 



(3.97) 



for a = 1, . . . , n. These are the GNLS equations |33 
from the AKS formalism (see also 
(loop algebra Q = Q <8>C[A, A -1 ] with Q 
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which have been derived in [Q entirely 
) associated with an afhne Lie algebraic structure 
sl(n + 1)). This construction reveals an affine Lie 
algebraic structure underlying the integrability of the SL(n + l,n) cKP hierarchy. 
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4 Darboux-Backlund Techniques of SL(p, q) cKP Hier- 
archies and Constrained Generalized Toda Lattices 



4.1 Emergence of the Toda structure from Two-Bose KP System 

The two-boson KP system defined by the Lax operator L = D + = D + a (D — b)~ x 

is the most basic constrained KP structure. It belongs to the SL(2, 1) cKP system. 

We start with the initial "free" Lax operator = D and perform a following transfor- 
mation: 



L« = (qWd^ 1 ) D (^D- 1 ^ 1 ) = D+U® (ln$W)"j D- 1 ($«») 1 (4.1) 

which we call a DB transformation. The construction below is a special application of 
properties listed in Appendix B including eq. ( |B.9| ). 

Successive application of DB transformations leads to the following recursive expressions: 

L (k+i) = f$(k) m (k)-^ L (k) (§(k) D -i$(k)-^ = D + Q(k+i) D -iy(k + i) 
$(k+i) = $ (fc) (ln$( fc ))"+ ($(*)) 2 #(*) , = (^(k)y 1 ( 43 ) 

Introduce now 

fc = In $ (fc) -» = e 0fe fc = 0, . . . (4.4) 

which allows us to rewrite ( f4.3| ) as a (ordinary one-dimensional) Toda lattice equation: 

d 2 (f) k = e **+ 1- ** - e^-^- 1 (4.5) 

Related objects ip n are: 

which satisfy due to eq. ( f4.5| ) the following form of the Toda lattice equation: 

with ip n = for n < 0. Comparing (4/?) with Jacobi's theorem (p. 7] ) we find the Wronskian 
representation for ^ n : 

i) n = In W n [$, 9$, . . . , d 11 - 1 ®] with $ = $ (0) (4.8) 

Correspondingly $( fc ) acquires the form: 

(fc) = W k+1 [*,d*,...,&*] 

W k [^,d^,...,d k - 1 ^} { ' ; 

We recognize at the right hand side of ( [4.7| ) a structure of the Cartan matrix for A n . Leznov 
considered such an equation with Wronskian solution (in two dimensions) in ||2] . 

Hence, the solutions of the (ordinary one-dimensional) Toda lattice equations, with 
boundary conditions ip n = for n < 0, reproduce the DB solutions of the ordinary two- 
boson KP hierarchy (|4.9|) upon taking into account that $ = $ ( - ^ = exp (0 O ) — exp C0i)- 
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Note that we can write W n = r n with the r-function r n satisfying Hirota's bilinear 
equation for the Toda lattice: 



r n d 2 T n - (dr n ) 2 = r n+ iT n _i 



(4.10) 



on basis of 

We will now find a linear discrete system (Toda spectral problem), which leads to the 
above structure. First define: 



(jj(n+l) 



T n+l T n-l 



S n = d(\n& n+ V)=d(\n^-) 

As result of ( [4.1 0|) and their definition S n , R n satisfy the Toda equations of motion: 

d S n = R n +i — Rn 
d R n = R n (S n — S n _i) 

These equations can also be obtained as a consistency of the spectral system 

d^ n = $ n+1 + S n V n 

A# n = ^ n+ i + S n ty n + Rn^n-i 



(4.11 
(4.12 



(4.13 
(4.14 



(4.15 
(4.16 



which defines a so-called Toda chain system. 

As we will show this purely discrete system contains information about the underlying 
continuous structure. This structure is being revealed when we realize that the lattice jump 
n —>■ n + 1 can be given a meaning of the DB transformation. We start by rewriting ( [4.15 
as follows: 

d$ n = V n + 1 + S n V n ~ V n+1 = ef s "de-f s ^ n (4.17 
or taking into account (|4.14|) as 

^ n+1 = R n ef s ^d (R n ef s »->-y 1 * n = ^^(n)^ = T(n)V n (4.18) 

where $(n) = R n e-f Sn ~ 1 and T(n) = $(n)cKI> _1 (n) plays a role of the DB transformation 
operator generating the lattice translation n — ► n + 1. 

The remaining of the Toda spectral equation ( 4.16 ) can be given a form (with ty(n) = 



A* n =(d + R n (d- Sn^y 1 ) ^ n = [d + R n eI s ^d- l e-I 



i i 



d + $(n)a _1 *(ri 



L(n)m, 



(4.19) 



of the Lax eigenvalue problem of the two-bose KP Lax system with generic two-bose KP 
Lax operator L = D + QD~ lj &. Hence the Toda lattice spectral problem has been shown 
equivalent to the continuous cKP system possessing the symmetry with respect to the DB 
transformations. 
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4.2 On the DB Transformations of the SL(r + n, n) cKP Lax Oper- 
ators 

We shall here consider behavior of the general class of constrained KP Lax operators from 
SL(r + n,n) cKP hierarchy under an arbitrary DB transformation L = xDx 1 L\D^ 1 X 1 
where \ is an eigenfunction of the Lax operator L as given in ( |3.26| ). The transformed Lax 
operator reads: 



§ 



L + + x(d x (x- 1 L + x) >1 D- 1 ) x' 1 



i=i 

n 



A 



i=i 



(4.20) 
(4.21) 
(4.22) 

xDx ^Lj X (4.23) 
(4.24) 



From the above discussion we know that all involved functions are (adjoint) eigenfunctions 
of L ( |3.26|) , i.e., they satisfy: 



_d_ 

dti 



f = Llf f = X ,*i 



d h. 
— *i = -L^j 

Otk 



(4.25) 



We are interested in the special case when x coincides with one of the original eigenfunctions 
of L, e.g. x — Then $i = and the DB transformation ( [4.20 ) preserves the form ( ^.26[ ) 
of the Lax operators involved, i.e., it becomes an awio-Backlund transformation. Applying 
the successive DB transformations in this case yields: 

L (*) = r (*-i) L (*-i)( T (*-i))- 1 = ( L (*)) + £$f>£r 1 *f> , T« = (4.26) 



(fc+i) 



$ ( fc+ i) = / T(fc) ^ $ ( fc) ^ $ (, 



$1' 



i=l 
-1 



fc = 0,1, 



(fc+1) 



spy 1 & (*m 



n 



(4.27) 
(4.28) 
(4.29) 



Using the first identity from ( |4.26|) , i.e., L (fc+1) T (fc) = T (fc) L (fe) , one can rewrite ( |4.27|) in the 
form: 



whereas: 



>(*) = T (k-i) T (k-2) . . . T (o) ^^(o)^$(°)^ 
= T ( k -l) T (k-2) . _ . T (o) $ (o) i = 2. 



(4.30) 
(4.31) 
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Finally, for the coefficient of the next-to-leading differential term in ( p. 26 ) w r -2 — r Res L^ = 
rd 2 \nr, we easily obtain from ( |4.21| ) (with X — ^1) its &-step DB-transformed expression: 



f \ 



(k) 
r-2 



— U 



(0) 



)=^ln^ y = ^ln($r i )...$f 



(0) 



(4.32) 



Remark. Let us particularly stress that the eigenfunctions we are working with are not 
Baker- Akhiezer eigenfunctions of L from ( |3.26j ). Imagine, namely that we started with 
L = L + + ipD~ x ip* , where ip, tp* are BA functions. Choosing x — V* would result according 
to (RL23D and ( ^24|) m L = D. Hence in such case we would be able to transform away the 
pseudo-differential part of the cKP Lax operator by a finite number of DB transformations. 



4.3 DB Transformation and Eigenfunctions of the Lax Operator 



After seeing DB transformation in action in the simple cases shown in the first two subsections 
of this chapter we are now ready to review the basic properties of DB transformation from 
the point of view of preserving the form of the Lax evolution equation. Related material can 
be found in e.g. |19|, |51|| . 

Lemma. For arbitrary pseudo- differential operator A we have the following identity |5^/ : 
X" 1 ) , = XD X ~ 1 (A) + X D-\- 1 - X d x (x~\(A) +X )) D^ 1 (4.33) 

the transformed Lax 



(xDx-'AxD- 1 



For L satisfying a general KP-KdV Lax equation d^L 
operator L = TLT~ X will satisfy: 



d k L 



TL r jT~ l + {d k T)T~ x , L 



(4.34) 



Let $ be an eigenfunction of L i.e. = dk&, which enters the DB transformation through 
T = $Z}$ _1 . One verifies easily that for A = L~ and x = & equation ( 4.33|) becomes 



■TLT-T- 1 

Correspondingly ( |4.34j ) takes the form: 

d k L 



TLlT- 1 + (dkT)T- 1 



L 



+ > 



(4.35) 



(4.36) 



and we have therefore established: 

Proposition. The DB transformation with an eigenfunction $ preserves the form of the 
Lax equation (\2.J\) i.e. the DB transformed Lax operator satisfies the same evolution equation 
as the original Lax operator. 

For the pseudo-differential operator A = D + a\D~ l + ... we have 



Res (xDx~ l A X D- l X ~ l ) = (In x)" + Res (A) 



(4.37) 



For the case of the Lax operator L with the r-function satisfying ResL = d 2 In r we get 
therefeore | JH| : 

Proposition. Under the DB transformation with an eigenfunction $ the r-function asso- 
ciated to the Lax operator L transforms according to r — > f = $r. 
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4.4 Exact Solutions of SL(p, q) cKP hierarchy via DB Transforma- 
tions 

Armed with the Wronskian identities from Appendix B, we can represent the fc-step DB 
transformation ( [4.30 ) — ( 4.32Q in terms of Wronskian determinants involving the coefficient 
functions of the "initial" Lax operator 

r— 2 n 

L (0) = jjr + £ u (0) D l + J2 $f (°) (4.38) 

«=0 i=l 

only. Indeed, using identity (B.9| ) and defining: 

fc = 1,2, ... (4.39) 

we arrive at the following general result: 

Proposition. The k-step DB-transformed eigenf unctions and the tau-function ( \(.30[ )— 
of the SL(r + n,n) cKP system ( 3. 71 ) for arbitrary initial ( \4 . 38j ) are given by: 

*f) = w k+l[ ^\x^-^} (440) 



wM$! 0) , x (1) ,...,x (fc - 1} ] 

^i[$S 0) ,X (1) ,-,X N 

iy fc [$i 0) , x w,..., x (^)] 



™ W^-iM 0) y (1) y (fe_1) $ (0) 1 

r W = H/ fe [$f ) , X (1) ) ---,X M ]r (0) (4-42) 
where t^\t^ are the r-functions of L(°\ , respectively, and \ is given by ([(.Sty . 

Example Consider SL(2, 1) cKP hierarchy with the Lax operator: L (0) = _D+$ (0) -D _1 \l/ (0) , 
which serves as a starting point of the successive DB transformations: 

= D + & k) D~ l ^ {k) = T {k ~ l) L {k ~ l) (r [k ~ l) y l (4.43) 

cjj(fc) = j.(fc-i)_^(fc-i)^(fe-i) (4.44) 

By iteration we find: 

*(*) _ T (fe-i) . . . T (i) T (o) /7/-(o)\k $ (o)\ _ ^fc+i[^ (0) ; X (1) ; --- ; X (fc) ] r44 ,x 

q(k) = (^(k-vy 1 ( 4-46 ) 

Comparing two alternative expressions involving the tau-function r^: 

T (k) = $ (fc) . . . $ (D $ (o) T (o) = Wifc+1 [$S 0) , X W, . . . , X (fc) ]r (0) (4-47) 
and ResL^ = $( fe )i]/( fe ) = <9 2 lnr( fc ) we arrive at: 

$(0)^(0) = -d 2 \nW k [^°\x {1 \---,X {k ' 1) } (4.48) 

w fc+1 [<&(°), x {1 \ • • • , x w ]w*-i[* (0) > - , x {k ~ 2) ] 



+ 



w fc 2 [$(°),x( 1 ),...,x (fe - 1) ] 
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which generalizes the structure of the Hirota equation Q4.10 ) found for the case $^°- ) = = 
in subsection 3.1. 

Example: Construction of the SL(3, 1) cKP Lax operator, i.e., r = 2,n = 1 in ( |3.26| ) i.e. 

L = D 2 + u + A(D - B)' 1 = D 2 + u + ^ZT 1 ^ (4.49) 
starting from the "free" = D 2 . This example is pertinent to the simplest nontrivial 



string two-matrix model [11 



From the basic formulas for successive DB transformations ( |4.27| )-( |4.26| ), applied to 
( pH| ), we have: 

L (k) = D 2 + u (k) + ^(k) D -ly(k) ( 45Q ) 

L {k) ^ L {k+i) = T (k) L {k) ^{^y 1 ^ rp(k) _ ^(fc) y) (foikfy 1 ( 45 -q 

u« = 2Res L l 2 = 2d 2 x In r« = 2d 2 x In (s^ 1 ) ■ ■ ■ <S>®] 



$ (fc) = A (fe) e /B< fc » = 



^•($(ti)^ 1)+2 ^ ln K" 2) --- $ 



(0) 



$(*-!) 
$(fc-2) 



(4.52) 
(4.53) 
(4.54) 



M (o) = ; v j / (o ) = . $( o) = /■ ^ c(A)e €(A,m) ? e(A>{t} ) = Ax + y A i t (4.55) 

Jr 2n ^2 

where $ < -°' ) in ( |4.55|) is an arbitrary eigenfunction of the "free" = D 2 (the contour T 
in the complex A-plane is chosen such that the generalized Laplace transform of c(A) is 
well-defined). 

As a corollary from the above proposition, we get in the case of ( f4.50| ) : 



= T (k-i) . . , T (o) ^2fc $ (0) 
T (k) = w [$(0)^(0)^ .. jfl £(*-D$(o) 





'$(°),fi£$< >,. 








$(°),<92$(o),.. 


^2(fc-l) $(0) - 



(4.56) 
(4.57) 



Substituting fl4.56|) , (|4.57|) into Q4.52j )- (|4.54j ) we obtain the following explicit solutions for 
the coefficient functions of ( f4.49| ) : 



a. In 



jy 


$(°),^$( ),. 


. .^2^-1)^(0)" 




$(°),«92$(°),. 


. .,^ (n_2) $(°)" 



B (n) 
A (n) 

In the more general case of SL(r + 1, 1) cKP Lax operator for arbitrary finite r : 

r-2 

L = D r + Y,uiD l + ®D- 1 ^ 

1=0 



w 


$ (0) ,d2^°),...,d2n$(0)" 


W 


$ (0) , ^$ (0) , . . . , c}| (n " 2) $ (0) 


(w 


^°),d 2 ^,...,d 2 x in ~ 1] ^ 


) 2 



(4.58) 
(4.59) 

(4.60) 



(4.61) 
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which defines the integrable hierarchy corresponding to the general string two-matrix model 
(cf. Jll], pi), the generalizations of (^56|) and ((157]) read: 



$ (fc) =T {k-l) . _ T (0) (gk-r^(0) 



w 



$ (0) ,<^$ (0) ,...,<9£' r $ (0) 



w 



$(°),a$(°),...,ai fe " 1)r $( ) 



1 (fc) 

r 



Res L~ = d 2 Jnr^ 
where is again given explicitly by ( [4.55| ). 



r « = W [$ (0) , <9;$ (0) , . . . , ^ 1 )- r $(°) 



(4.62) 
(4.63) 



4.5 Relation to the Constrained Generalized Toda Lattices 



Here we shall establish the equivalence between the set of successive DB transformations of 
the SL{r + 1, 1) cKP system (WM ■ 



L (k+l) _ T (k) L (k) ^ T (k)^ 1 ^ T {k) _ §(k) D $ (fc)- 1 

L (0) = £,r + g u (0) D J + $ (0) D -1^(0) 



(4.64) 
(4.65) 



and the equations of motion of a constrained generalized Toda lattice system, underlying 
the two-matrix string model, which contains, in particular, the two-dimensional Toda lattice 
equations. 

For simplicity we shall illustrate the above property on the simplest nontrivial case of 
SX(3, 1) cKP hierarchy flQgp . We note that eqs. fl£33)- fl£54p (or (|g|)-(g6g)) can be cast 
in the following recurrence form: 



u 



2d T B {n) 



with "initial" conditions (cf. ( |4.55| )) : 

A (o) = B (o) 



u 



(0) 







B w = SL In $ 



(4.66) 
(4.67) 

(4.68) 



(4.69) 



where $ is so far an arbitrary function. Now, we can view ( |4.66| )— (]4.68|) as a system of lattice 
equations for the dynamical variables A^ n \ B^ n \ associated with each lattice site n and 
subject to the boundary conditions: 



U 



(n) 







n < 



(4.70) 



Taking ( |4.69| ) as initial data, one can solve the lattice system (|4.66|) - (|4.68|) step by step (for 



n = l,2, 



.) and the solution has precisely the form of (|4.58|) -( |4.60|) . 
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4.6 Darboux-Backlund Transformation and the Dressing Chain. 

As a small digression we will here apply Darboux-Backlund (DB) transformation to the KdV 
hierarchy emphasizing similarity with the approach developed in the KP setting. Recall from 
fl£5q)-(l£55D that for T = QDQ' 1 we have 

TD 2 T- X = D 2 + 2<92ln$ + $($- 1 $") , /r 1 $- 1 (4.71) 

The inverse scattering problem corresponding to the KdV equation d t u + 3udu + ^d 3 u = 
is defined in terms of the differential operator L = D 2 + u, which transforms as: 

L X = T [D 2 + u) T- 1 = D 2 + 2d 2 x In $ + $ (VV)' D' 1 ^' 1 

+ u + $u' 'D- 1 ®- 1 (4.72) 

In order for L\ to be a differential operator we have to demand that the pseudo-differential 
part in ( |4.72j ) vanishes: 

+ W = (4.73) 



which translates into relation between u and $: 

^ = -^-A = -/o 2 -/o-Ao (4.74) 
where Ao is an integration constant and fo = (ln$)'. In this notation: 

L = D 2 + u = D 2 - f 2 - f Q - A = (D + f ) (D - /„) - A (4.75) 
Since T = D — f and = —D — f we can factorize the Lax operator L according to 

L = -r t T + A (4.76) 
In this notation the DB transformation on L: 

L x = TLT- 1 = -TT^-X = D 2 - f 2 + ft - A = D 2 + Ul (4.77) 
amounts to reversing the T, T' operators in expression for L. Moreover we find that 

ui = -/ 2 + ft - Ao = - A 2 -f[~Xi = u + 2(ln$)" (4.78) 

where we introduced new variables f\ = (ln$i)' and Ai, which allow us to rewrite L\ as 
L x = -T{ Ti + Ai with Ti = D- f\. Clearly we have Li$ a = -Ai$i. 

Successive use of the DB transformations |55| leads to the Lax operators L n = D 2 + u n 



with 

U n = -ftl + fn-l ~ An-1 = "/* " & " A n = U + 2(ln $ • • • S^)" (4-79) 

with a string of eigenvalue problems: 

L&i = -\i<$>i or (d 2 + m + A^ $i = $ = $ , m = u (4.80) 
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for i = 0, 1, Let *ffx, . . . , be solutions of the related eigenvalue problem: 

= or (<9 2 + u + Ai) ^ = = $ , i = 0, 1, . . . (4.81) 



We have a following theorem shown by Crum [23], which establishes relation between <3>j and 

eigenfunctions. 
Proposition. The function 



$ = 



W w+ i[^q,^i,...,^ w ] 



satisfies the differential equation: 

(d 2 



U, 



u n + A n ] <i.\, = 

J2 



M + 2«9 2 In W n [* , *n-l] 



(4.82) 



(4.83) 
(4.84) 



Proof. To prove it let us notice that L n = (T n _i • • • To)Lo(T n _i • • • To) -1 , with Tj = 
^fid^fj' . From L n $„ = — A n $ n we find that $ n = T„_i • • -To^. By use of induction 
and ( B.9|) one completes the proof. Equation fl4.84j ) follows automatically from (|4.82| ) and 

(EZD- 

Note a clear analogy of the dressing chain construction with the Susy QM j|7|]. Define 
namely, T n — D — f„ with = — (D + f n ) and introduce H± through: 



T n Ti = -(D 2 -f 2 n + f' n ) = H_-\ n ■ TlT n = -{D 2 -fl-f' n ) 
These relations can be cast into the QM Susy algebra: 

{Q-,Q + } = H-\ n ; [Q±,H] = 
in terms of the 2x2 matrices: 



Ha- — A r 




T 



Q- 



Tt 




H 



H 

H, 



(4.85) 



(4.86) 



(4.87) 



The eigenvalue problem HV n = \ n V n with V n = (V_, V + ) T takes a familiar form (.D 2 — / 2 J 
f' n (Jz)V = and solutions are given by 



V ± = exp (± J / n ) = ($ f 



(4.88) 



as verified by rewriting the eigenvalue problem as T^TV + = and TT^ VI = 0. Hence the 
Darboux techniques prove useful in constructing exact solution of the Schrodinger problems 
of the supersymmetric quantum mechanics. 
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4.7 Connection to Grassmannian Manifolds and n-Soliton Solu- 
tion for the KP Hierarchy 



Let {^i, . . . ,ip n } be a basis of solutions of the n-th order equation Dip = 0, where L = 
(D + v n ) (D + v n -i) . . . (D + vx). If Wk denotes the Wronskian determinant of {^i, . . . , ipk} 
one can then show that |60|, 1381 : 



Vi = d{hJ^j W = l (4.89) 

This allows to establish that the space of differential operators is parametrized by the Grass- 
mannian manifold (see e.g. 46|). Start namely with the given differential operator 



L n = D n + u\D n ~ 1 + • • - + u n and determine the kernel of L n given by n-dimensional subspace 
of some Hilbert space of functions Ti, spanned, let say, by {ipi, . . . ,ip n }- This establishes 
the connection one way. On the other hand let {ipi, . . . ,ipn} be a basis of one point fl of 
a Grassmannian manifold Gr^ n \ Define the differential equation as L n (Q)f = Wk(f)/Wk- 
From (|B.11|) this associates the differential operator 



L n = D n + u x D n - x + • ■ ■ + u n = (D + v n ) (D + v n ^) ■ ■ ■ (D + Vl ) (4.90) 

given by a Miura correspondence to a given point of the Grassmannian. 

Let us now comment on connection to n-soliton solution for the KP hierarchy. Assume 
that the above functions ipi i = 1, . . . , n have the property d m ipi = d m ipi for arbitrary m > 1 
(d m = gf-), in other words ipi are eigenfunctions of = D. We introduce L = LnDL^ 1 , 
where L n is defined in terms of {^i, . . . , ip n } as in ( |4.89| ) and QL9Q|) . It is known |2^, 



that such a Lax operator satisfies a generalized Lax equation d m L = [L™, L\. 

Using ( [B.9| ) we can rewrite the above Lax operator as a result of successive DB transfor- 
mations applied on D: 

L = L n DL~ l = T n T n _! • • • Ti DT{ 1 ■ ■ ■ T~l x T~ x (4.91) 

where Tj are given in terms of Wronskians as in ( |B.10| ). It follows that L can be cast in 
a form of the Lax operator belonging to subclass of the SL(1 + n,n) cKP hierarchy and 
having the form as in ( |3.26| ) with r = 1. Using the formalism developed in this paper 
one can prove by induction that the corresponding r-function of L takes a Wronskian form 
Tn = W n [ipi, . . . ,ip n ] reproducing n-soliton solution to the KP equation derived in fl33| . In 
fact, choosing ^ = exp (x^fcQ^) + ex P (S^fc/^f) allows to rewrite r n in the conventional 
form of the n-soliton solution to the KP equation ^ . 
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5 Two-Matrix Model as a SL(r + 1, l)-cKP Hierarchy 

5.1 Two-Matrix model, Orthogonal Polynomials Technique 

We shall consider the two-matrix model with partition function : 

„ r pi P2 ~\ 

Z N [t, t,g} = J dM 1 dM 2 exp j £ t r Tr M{ + ^ t s Tr M s 2 + g Tr M X M 2 \ (5.1) 



where M 12 are Hermitian N x N matrices, and the orders of the matrix "potentials" pi )2 
may be finite or infinite. As in ref.|0| and |57|, |J we will use the method of generalized 
orthogonal polynomials [|T7| to evaluate partition function (5TT). After angular integration 
in (|5.1|) we obtain 

I r N ~ i N ~ 1 

Zn% i,9] = j7j j]l d\d\ exp £ (y(\) + V(\i) + g\ty A(A t ) A(A<) (5.2) 

where V(A) = E4A fc , V(X) = E^A 5 . A(Aj), A(\) are standard Van der Monde determi- 
nants. As in one-matrix model one can deal with them using orthogonal polynomials. Here 
we will work with two families of the orthogonal polynomials: 

Pn(Ai) = A^ + 0(Ar 1 ) ; J R m (A 2 ) = A^ + 0(A™- 1 ) , n,m = 0,1,... (5.3) 
which enter the orthogonal relation: 

J dX x d\ 2 exp (V(Ai) + K(A 2 ) + <?AiA 2 ) P n (A!)P m (A 2 ) = /i m 5 nm (5.4) 

Following approach to one-matrix model we write down the recursion relations for the poly- 
nomials in (|5.3|) . 

n+l n+l 

A 1 P n (A 1 ) = £ Q^P(Ai) ; A 2 P n (A 2 ) = £ £(A 2 )Qin (5-5) 

From the definition of orthogonal polynomials it follows that 

Qn,n+i = 1 ; Qn,i = l>n + 2 and Qn+i,n = 1 ! <5i,n = l>n + 2 (5.6) 

The orthogonal polynomials approach leads to expression for the partition function in terms 
of h n : 

Zn — hohi ■ ■ ■ /ijv-2^jv-i (5-7) 
From the orthogonal relation ( |5.4|) and definitions (|5.6|) we obtain: 

= ; -ar = W- (5 - 8) 



36 



as well as 



dP n 

dt r 

dPn 

dt s 

where we introduced 



n-l 



£<&fl(Ai) ; 

1=0 

n-l 



1=0 



OPn 

dt r 

dPn 

dL 



n-l 

1=0 

n— 1 
Z=0 



Qnm = [HQH J 

Define now a wave function 



with H nm — h n 5 nm 



^ n (t,t,X)=P n (X 1 )e v ^ 



(5.9) 
(5.10) 

(5.11) 
(5.12) 



From definition of \l/ n and ( |5.5|) , ( 5.9|) and ( |5.10|) we obtain the following eigenvalue problem: 



(5.13) 
(5.14) 

(5.15) 



where \I> is a semi-infinite column (. . . , ^n+ij • • •) an d Q and Q are semi-infinite matrices, 
i.e., with indices running from to oo, Furthermore we adhere to the following notation: the 
subscripts — /+ denote lower/upper triangular parts of the matrix, whereas (+)/(—) denote 
upper /lower triangular plus diagonal parts. 
In addition we also have a relation: 



saT 



(5.16) 



which can be derived using identity: 

d 



= jdXxdXi— e^(y(X 1 )+V(X 2 )+gX 1 X 2 )P n (X 1 )P m (X 2 ) i = l,2 (5.17) 
5.2 Two-Matrix Model as a Discrete Linear System 



The compatibility of the eigenvalue problem ( |5.13| )-( |5TT6| ) [ p~5| , |TTf gives rise to a discrete 
linear system which we shall call the constrained generalized Toda lattice hierarchy: 



d_ 

dt r 



Q — Q(+), Q 



d_ 

dt r 



Q — Q(+) > Q 



-^rQ = [Q, Qi 

Ot s 



^rQ=[Q,Q- 

Ot s 

-g \Q, Q = 



r = 

8 = 1,.-.,P2 



(5.18) 

(5.19) 
(5.20) 
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As will be shown later the lattice system ( }4.66| )-( f4.68| ) can be identified with the t\ = x 
evolution equations of the above discrete linear system and therefore provides a direct link 
between matrix models and the integrable hierarchies PS|| . 

In what follows it will be convenient to define an explicit parametrization of matrices Q 
and Q. We choose the following parametrization which is consistent with G5.6|) : 



Qnn = ao{n) , Q n ,n+X = 1 , Qn,n-k = Ofc(n) k=l,...,p 2 ~l 

Qnm — for m — n > 2 , n — m > p 2 (5.21) 

Qnn = bo(n) , Qn,n-1 = Rn , Qn,n+k = ( n )R n +l ' ' ' R-n+k k = l,...,Pl~ 1 

Qnm — for n — m > 2 , m — n > pi (5.22) 

In most examples used in this section we work with the number p 2 = 3 in (|5.21| ), whereas 
the number pi in ( |5.22j ) remains finite or oo []. 

Let us consider for the moment the first evolution parameters ti,ti as coordinates of a 
two-dimensional space, i.e., t\ = x and ti = so all modes ak{n),bk(n) and _R„ depend on 
(x, 2/;i 2) • • • , t Pl ;t 2 , . . . , t P2 ) . 

The second lattice equation of motion ( |5.19| ) for s = 1, using parametrization ( |5.22| ), 
gives : 

d*R n = R n (b (n)-b (n-l)) , <3 5 &o(n) = &i(n) - b x {n - 1) (5.23) 
flfef W p ) = ^ + i(n)-6 fc+ i(n-l) ) fc > 2 (524) 

\-Rn+l • • • Rn+k J Rn+1 ■ ■ ■ Rn+k 

Similarly, the second lattice equation of motion (|5.18|) for r = 1 gives : 

d y b Q (n) = R n+1 - R n , d y b k (n) = _R n+ i6fc-i(™ + 1) - i? n+fc 6 fc _i(n) (5.25) 

for k > 1. From the above equations one can express all b k {n ± £) , k > 2 and R n ±i - 
arbitrary integer) as functionals of b (n), bi(n) at a /ixed lattice site n. 

Furthermore, the lattice equations of motion ( |5.18| )-( [5.19| ) for r = 1, s = 1 read explicitly: 

d s a (n) = R n+l - R n , <9 £ a fc (n) = R n _ k+ ia k ^i{n) - R^^n - 1) (5.26) 

(with A; > 1) and 

a / \ / i i \ / \ o / afcW \ a fc+ i(n+ 1) - a fc+ i(n) 

c\a (n) = ai(n + 1) - ^(n) , d„ = (5.27) 

\JrC n . . . JrC n -k+l J -tt n ■ ■ ■ -ttn-k+1 

with k > 1. Following (pT23ft , ( ^25| ), ( gT26|) we obtain the "duality" relations: 

d y bi(n) = dxR n+ i , <9 £ a (n) = d y b (n) , d^a^n) = d y R n (5.28) 
From the above one gets the two-dimensional Toda lattice equation: 

d y \nR n = a (n) - a (n - 1) — > <9zc\ lni? n = - 2/? n + (5.29) 

3 Both numbers pi.2 indicating the number of non-zero diagonals, outside the main one, of the matrices 
Q and Q are related with the p olynomial orders of the corresponding string two-matrix model potentials, 
whereas the constant g in ( 5.20 ) denotes the coupling parameter between the two random matrices. 
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Eqs.( |5.26[) -( |5.29[) allow to express all akin + £) , k > 1 and R n ±e as functionals of a (n) 
and R n (or ai(n) instead) at a fixed lattice site n. Furthermore, due to eqs. fl5.26|) and ( |5.25|) 
for k — 1, all matrix elements of Q and Q are functionals of bo(n),bi(n) at a fixed lattice 
site n. Alternatively, due to ( |5.28| ) we can consider a (n) and Rn+i as independent functions 
instead of bo(n), bi{n). 

Let us also add the explicit expressions for the flow eqs. for b (n), bi(n), R n+ i resulting 
from fU8l) and (PD : 











dt 
d_ 

dt, 



■bo(n) = d*(Q r 



nn Qf 

d_ 

)nn ' dt r 







r b (n) = d s (Q s ) , -^(n) = dJR n+1 (Q s ) ] , w R^i = d s (Q s ) (5.30) 

f_ \ J nn (jt „ L V /n,n+V Qf g \ / n+X,n 

hin) = d,[R n+1 (Q r ) n , n+ i] > ^Rn+i = d £ (Q r ) n+hn (5.31) 



5.3 String Equation 

Note the presence of the non-evolution constraint equation ( |5.20|) , which is called a string 
equation. The lattice equations for the matrix elements afc(n) of Q (the first eqs. ( p. 18] ) and 
Q5.19Q ) can be solved explicitly as functionals of the matrix elements of Q : 



P2-1 

Q(-) = E a *QU 

s=0 



O'.s- 



(8 + 1) 



(5.32) 



Equation (|5.32|) follows simply from (|5.17| ) with i = 2. 

Note also that there is a complete duality under pi < — > p2 when the order pi of the first 



matrix potential in ( |5.1|) is also finite. For instance due to this duality we can obtain the 
analog of eq. ( p.32| ): 



pi-i 

Q(+) = E PrQ\+) 
r=0 



Pr 



(r + l) 



(5.33) 



P2 , t s < — > t r , Q(-) < — ► Q(+)- Equation (|5.33| ) follows simply 



by interchanging pi <- 
from ( 5.17 ) with i — 1. 

The first equations of (|5.26| ) and (|5.29|) , imply the following two additional constraints 
on the independent functions b (n) and R n+ i (or &i(n)) 



d y b {n) = d x I a sQ S r 



dyRn+l = 9 X a sQn+l,n 



(5.34) 



s=0 



s=l 



In fact, eqs. (|5.34 ) are nothing but the component form of the string equation ( 5.20 ). Indeed 
upon using ( |5.32|) and ( |5.18|) , ( |5.19| ), eq. ( |5.20| ) can be rewritten in the following form: 



P2_1 \ - 1 

d y -J2 Usd/dt s Q = — 1 
s=i ) 9 



(5.35) 
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or equivalently 



9 v - X a sd/dt s \b (n] 

s=l 

P2-1 

dy ~ X a sd/dt s ] = 

s=l 



1 





P2— 1 



^ - ^ a s 0/0* s 6*(n) = (5.36) 



8=1 



for fc > 1. Now, inserting ( |5.30| ) into ( 5.34 ) we find that the latter two equations precisely 
coincide with the nn and n+1, n component of the string eq. ( |5.36|) . Hence the string equation 
would amount to identifying the flow d/dy = d/dt\ with the flow 







Of 



P2-1 «=1 



(5.37) 



if not for the constant — 1/g on the right hand side of (|5.36|) . For this reason it is more 
convenient to make a change of variables and correspondingly introduce another matrix Q: 



QT-) 1 = E *sQU 



P2-1 



s=Q 



We now generalize (|5.38|) to: 



QL) = £ IsaQ?-) s=i,...,p 2 



(5.38) 



(5.39) 



cr=0 



and parametrize Q as in ( |5.22| )) with matrix elements R n ,bk(n) in place of R n ,bk(n). Coef- 
ficients jso- are simply expressed through a s = 7 P2 -i,s (7oo = 1) as: 



7ss = (Til)' , 7s,s-l = - (7ll) S 2 721 , 1s,s-2 = (in] 



8(8-3) , - 2 . 8 

g (72i) + -7n73i 



(5.40) 



7n = {a P 2-i) p ' 



721 = 



a P2 -2 



731 



Q!p 2 _3 



P2-4 



a 



P2-2 



P2 - 1 

and for the Q-matrix elements we obtain: 

721 



2P2-5 



(5.41) 



i? n = 7n-R n , fe (n) = 7n & o(n) + 



2 7 



61 (n) = 7iA( n ) + 



731 



721 



11 



37ii V 2 7 



(5.42) 



11 



etc.. 



We can now rewrite the string equation in the "hatted" variables from eq. (|5.42|) . First 
we note that the explicit calculation based on ( p.41[ ) gives 



d s _ <9 s + 1 
7*,--i = — ls,s and 7 SjS _ 2 = — 7 SjS _i 



(5.43) 
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Especially, we find 



and 



Hence we easily get 



d 



dt 



-721 







P2-1 



731 



2 2 



dt 



P2-1 



721 



37 



11 



27 



11 



<9 



5 



= 



i = 1,2 



\<9tl Sipa-i, 

and the same result for R n . This allows us to postulate a general statement: 

( d d \ a 

Q = 



9ii <9t„ 2 _i 



(5.44) 



(5.45) 



(5.46) 



(5.47) 



Accordingly in the new "hatted" framework the string equation boils down to identification 
of the flow d/dti with d/di P2 ^i flow. 

Consider now an arbitrary matrix X(Q, Q). From equations fl5.19|) , ( |5.18| ) and Q5.32D we 
find 



dX(Q,Q) 



dt 



[X,Q. 



P2-1 



Q W ,X + [X,Q\ 



dx 



+ [X,Q[ 



(5.48 



Comparing with ( |5.47| ) we find that Q , Q = 0. It is therefore plausible to supplement 
( p.38| ) by condition 

g^y 1 = /+ — q^ 1 = q (5.49) 

with {I+) nm = bn+x,m- The last equality follows from ( |5.32| ). Let us also introduce the new 
subset of evolution parameters {4} through the equations: 



s = l,...,p 2 



(5.50) 



which for s = p 2 — l reproduces (|5.37|) . Now, taking into account (|5.38Q , ( |5.50|) and condition 
( |5.49j ), all constrained Toda lattice eqs.( |5.18| )-( p.20| ) can be re-expressed as a single set of 
the flow equations for the matrix Q : 



_d_ 

dt. 



Q= Q,Qt 



S = 1,...,P2,2(P2-1),3(P2-1),---,P1(P2- 1] 



with the identification t r = t r ( P2 _i), The associated linear problem reads: 





dt 



■A 



A, 



(5.51) 



(5.52) 



From the fact that Q has the same form as Q up to the "hatting" we can deduce that Q 
satisfies the Zakharov-Shabat equation: 



ot r at* l 







(5.53) 
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Recalling associated linear problem ( |5.52| ) we see that 

i> = o 



translates into: 



3_ d_ 

dL ' dt , 



(5.54) 



establishing the commutativity of flows. In fact an explicit calculation for the simplest 
nontrivial case of p 2 = 3 gives indeed ( |5.54| ) for r, s — 1, 2, 3 ! 

Recall now the equation ( |5.33| ). Together with the string equation ( p. 2 Op it imposes a 
following condition on Q: 



'pi— l 



5 



'pi-i 



d 



-i ^ 



r=l 



r=l 



dt r 



9 



(5.55) 



Due to identification between Q and Q modes in (|5.49 ) this extra condition translates into 
a condition for Q P2 ~ l . 

We obtain from ( |5.51P ,( 5".52|) the following continuum Lax problem at a fixed lattice site 

n: 



4-L(n) 



L(n)ip n 



L(n) , £ s (n) 



ip n = -C s (n)ijj n (5.56) 
, s = l,...,p 2 ,2(p 2 -l),3(p 2 -l),...p l (j> 2 -l) (5.57) 



L{n) = -D- 1 R n+1 + OST 1 + £ 



Pa-l ('_Ti fc nP 2_1 
1 . v ' ^n,n-k 



k—1 1L n 



R n . ■ ■ R 



n— fc+1 



Ac — <9rc In (-Rn • • • Rn-k+2 ) ) X • • ■ 



. . . X 



£ s (n) 



k=l Rn ■ ■ ■ R n—k+Y 



\D X — d x In (h n . . . Rn-k+2) ) • • • D x 



(5.58) 
(5.59) 



where all coefficients are expressed in terms of R n+ i, b (n), . . . , & pi _ 2 (n) at a fixed site n 
through the t% = x lattice equations of motion (|5.51|) . Up to a gauge transformation and 
conjugation, the explicit form of L(n) reads: 



e fbo(n) e -Jio(n) = jjn-x + (p2 _ x) h x {n)Dl*-* + ■■■ + R n+1 (D x - S (n 



-i 



(5.60) 



Eqs . ( |5 . 5 6| ) - ( |5 . 5 9| ) are the continuum analogs of the constrained Toda lattice Lax equations 
Q5.13] )-( 5.i6| ) and (|5.18|) -( 5.20|) without taking any continuum (double-scaling) limit. Let us 
particularly stress that they explicitly incorporate the whole information from the matrix- 
model string equation (|5.20| ) through (|5.32|) and ( |5.38| )- (|5.42| ) which were used in their 
derivation. 

Consider a case p 2 = 3. The Lax operator L(n) from ( |5.58|) is given by: 



L(n) = -IT 1 R n+1 + 2h(n) + b 2 n {n) - d x b {n) - 2b (n)D + D z 



(5.61) 
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The corresponding lattice equations then read (we write down explicitly only the i\ = x and 
t 2 evolution equations ) : 



d x lnR n+ i = b (n + 1) - b Q (n) , h(n) - bi(n - 1) = d x b (n) 
R n+1 -R n = d x (b 2 {n) + &i(n) + h{n - 1) 



d * r « a >• 

TT^-Rn+i = <9ic d x R n+ i + 2b (n)R n+ i 

_6 (n) = a, 26i(n) + 6g(n) - <9x6 (n) 
at 2 L 



9 

-5-61(71) = <9 x i? n+ i 



(5.62) 
(5.63) 

(5.64) 
(5.65) 



Now, we observe that the system of Darboux-Backlund equations for SL(3, 1) cKP hierarchy 
( |4.66j )-( |4.68j ) exactly coincides upon identification: 



B^ = b (n-1) , «W = 26i(ti-1) 



= R r 



(5.66) 



with the x = t\ constrained Toda lattice evolution equations ( |5.62| )- fl5.63p . Also, the higher 
Toda lattice evolution parameters can be identified with the following subset of evolution 
parameters of the SL(p 2 , 1) KP-KdV hierarchy flOH] ) JTT], : 



- _ ,KP-KdV 1 



(5.67) 



the second identification resulting from ( [5.38 



In particular, excluding 60(71) = B^ n+1 ^ and 61 1 



n 



in ( ggg ) using (gB2|)-(ggg), 



we obtain the two-dimensional Toda lattice equation for A^ n ^ = i? n : 

4% In A (n) = A {n+1) - 2A {n) + A {n ~ l) 

5.4 Partition Function of the Two-Matrix String Model 



(5.68) 



The partition function of the two-matrix string model is simply expressed in terms of 
the Q matrix element b\(N — 1) at the Toda lattice site N — 1, where N indicates the size 
(N x N) of the pertinent random matrices: d 2 x InZ^ = 61 (iV — 1) (cf. |L5|, Pfl). Thus, using 
( |4.42| ) and ( |4.55| ) together with ( |5.67| ), and accounting for the relations (|5.42| )-( |5.41|) , we 
obtain the following exact solution at finite N for the two-matrix model partition function: 



Z N 
$(°) 



W. 



N 



(0) 



exp — N 



721 



27 



11 



- C (A)e«(MM}) 



pi 



(5.69) 
(5.70) 



where W N [...]= det (dt 1 ^ 1 ¥°A 

* V p 2-i / l<i,j<N 



e (a, {*,*}) = £ A ^ + E 

s=l r =2 

with x = ti. The 7-coefficients are defined in 



( |5.41| ). The "density" function c(A) in Q5.70| ) is determined from matching the expression for 



43 



$(°): ^lnl>( ) = 6o(0) = 7n&o(°) + 72i/27n (cf. ($Mj) , with the expression 
for 6q(0) in the orthogonal-polynomial formalism fll5| : 



I PI 

d\id\2 exp I ^2 ^l^r + E + S , AiA 2 



(5.71) 




P2 

E 

s=l 



One can easily check that Wn from (|5.69|) satisfies the two-dimensional Toda lattice equa- 
tion d 2 IiiWn I dt\dt P2 -i = Wn+xWn-x/Wj^. Furthermore Z^ii^ provides a solution 
to the two-dimensional Toda lattice equation d 2 In Z^jdt\dt\ = Zjv+i^jv-i/ underlying 
the two-matrix model. This shows relevance of methods based on the Darboux-Backlund 
transformations for obtaining partition functions of the multi-matrix models. 
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A Schur Polynomials 

Definition. The Schur polynomials p(xi,x 2 , ■ ■ .) are polynomials of the multi- variable x = 
(xi, X2, ■ ■ •) defined by the generating function: 

(oo \ oo 

Y J x k z k \=Y J Pk{x)z k (A.l) 
k=l I k=0 

and Pk(x) — for k < 1. 
In particular 

p (x) = 1 , px(x) = xi , p 2 (x) = xj/2 + x 2 , p 3 {x) = xf/6 + x x x 2 + x 3 (A. 2) 

and generally p n (x) = x™/n\ + . . . + x n . Also (d/dx m )p n (x) = p n - m {x). 

We need the following recurrence relation. 
Lemma. 

m 

mp m (x) = X kPm-k (%) (A. 3) 



k=l 



where x = (xi, . . .) 

Proof. From definition (|A.1|) we get 



\fc=l K ) k=0 



Multiplying on both sides by J2Ho x i z i we obtain 



oo / m 



[ x kPm-k (x))z m = f^xiziew(it^ z l = z 4- ex P ( £ ~T zk 

m=0 \k=l / 1=0 \k=l K ) aZ \k=l K / 



d 00 00 
Z ~T E Prn(i)z m = Y, mp m (x)z m (A.5) 

® Z m=0 m=0 



the final result follows now by comparing coefficients of z Tl 
Especially taking x k = —d/dtk in ( |A.3| ) we get 



, ~ _ ™ a Pm _ fc (-0) _ ^1 d Pm ^ {-§) Q 



(A.6) 



which reproduces (|2.21|) in the text when applied on In r after an extra x-derivative is taken 
on both sides. 



B Wronskian Preliminaries 

We list here three basic properties of the Wronskian determinants. 
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1) The derivative T>' of a determinant T> of order n, whose entries are different iable functions, 
can be written as: 

V = V (1) +V {2) + ... + V (n) (B.l) 
where T>m is obtained from D by differentiating the entries in the i-th row. 

2) Jacobi expansion theorem: 

W k (/) W k .x = W k W' k _ x (/) - W' k W k . x (/) (B.2) 



or 



where the Wronskians are 



W k = W fc [Vi, . . . , Vfc] = det { (B.4) 



l<i,j<* 



and 

Pn_r(/)EE W k [^...^ k - U f]. (B.5) 

The proof (see also [fj, ||) of ( [B. 2| ) uses the fact that the left hand side of (|B.3| ) can be written 
as a linear combination Y^=% a pf^- Expression for the coefficients a p (up to multiplication 
with a common function) can be obtained by e.g. the Cramer's rule due to the fact that 
J2p=\ QpW — for i = 1, . . . , k as can be verified directly from ([B.3Q . It is now easy to see 
that the left hand side must be proportional to W k (/). The proof follows now by establishing 
that the terms with on both sides agree. 

Take now a special class of Wronskians 

W n = W n ty, V', • • • , d n ~ V] = det (^^) ■ (B.6) 



l<i ,j<n 



Hence, from (P-2|) we get : 

W n W'; - W' n W' n = W n W' n _ x («9» - W' n W n . x (<9» = W n ^W n+1 (B.7) 
which can be rewritten as 

&\*W n = Wn+ ^ (B.8) 



n 



3) Iterative composition of Wronskians: 



T fe T fe _! •••T 1 (/) = (B.9) 

vV k 

where 
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The proof of (EO) follows by simple iteration of ( |B.3| ) (see also the standard references on 
this subject [^3|, [3£| §]). For future use let us rewrite ( P-9| ) as: 



(a + u fc )(9 + i; fc -i)---(9 + T;i) / 



Wkif) 



Wj-i 



Vj = <9 r In 



(B.11) 



or 



(£> + Ufc) (-D + ...(£> + Vl ) = _ 



^1 ^2 



1 

D 

D k 



(B.12) 
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